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We have technically improved the non-perturbative renormalization method, proposed by Martinelli et al., by 
using quark momentum sources and sinks. Composite two-fermion operators up to three derivatives have been 
measured for Wilson fermions and Sheikholeslami-Wohlert improved fermions in the quenched approximation. 
The calculations are performed in the Landau gauge on 16 3 32 lattices at ft = 6.0 for 3 n values in each case. The 
improved sources greatly decrease the statistical noise. We extract and discuss here renormalization factors for 
local operators and moments of the structure functions for Wilson fermions. 



1. INTRODUCTION 

In this paper we want to discuss some gen- 
eral and technical aspects of calculating non- 
perturbatively renormalization factors of bilin- 
ear quark operators by imposing renormalization 
conditions on off-shell quark Green's functions on 
the lattice §. 

This method offers the possibility of comput- 
ing non-perturbative contributions and all orders 
of QCD perturbation theory in a simple way in 
contrast to the enormous efforts to get only lower 
orders with perturbative methods. The disadvan- 
tages of this numerical renormalization are the 
need of computer power and at present the sys- 
tematic uncertainties due to a quenched approxi- 
mation (that excludes fermion loops), discretiza- 
tion errors Q and gauge fixing. 

We look here at operators || O ~ qOq with 
quark fields q and O constructed out of Dirac op- 
erators 7 M and covariant lattice derivatives D v . 

2. METHOD 

The renormalization method we study is a 
M OM-scheme using the lattice as regulator with 

*Talk given by H. Oelrich at Lat97, Edinburgh, U.K. 



the lattice spacing a as cut-off. Lattice operators 
0(a) and quark fields q(a) are renormalized be- 
ing then only a function of the scale parameter 
via: 

O*0i) = Z ({an) 2 ,g(a))0(a), 
q R {n) = zV*((atf,g(a))q(a) 

(g{a) is the bare lattice coupling.) 

To determine the renormalization constants Zq 
one imposes conditions on amputated forward 
quark vertex functions in momentum space with 
external four-momenta p: 

(q{p)\0^)\q{p)) B amp = (q{p)\0{a)\q{p)Y a Z e p 

With a projector method, where Tr[To x ...] is 
taken, Tr being a colour x spin trace, these pre- 
scriptions give: 

z o Z q 1 Tr\r {q{jp)\0{a)\q{p)) 

= Tr[T (q(p)\0(a)\q(p))lZ e p }\ p2= ^ M 

As an optimal choice for the projector, we use, 
written in a general form: 

r ~(g(p)|0(o)|g(p))S^ 

There are two definitions for Z q : Usually 
one projects onto the energy-momentum part 
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sin(p M a)7^/a (continuum form: of the inverse 
propagator S 1-1 and gets: 

(The normalization is chosen to give Z q = 1 in 
the free case.) Another definition for Z q can be 
derived from eq. (1), setting O equal to the 
conserved vector current J M and using Zj = 1. 
But this is only compatible to the first deter- 
mination of Z q , if the component of J M trans- 
verse to p is taken. (This can be proven by a 
Ward Identity J^(p,p) = —id/dp fl S^ 1 (p) assum- 
ing S^ 1 (p) = iZ q (p 2 ) p"+ B(p 2 )l, which is true in 
the continuum and approximately true on the lat- 
tice for not too large p 2 a 2 .) 

3. NUMERICAL IMPLEMENTATION 

In a first step Monte Carlo gluon field config- 
urations U have to be generated and numerically 
fixed to the Landau gauge jl| . The lattice matrix 
elements in eq. (1), non-amputated, can then be 
calculated in those background fields: 

(q(p)\0(a)\q(p)) 

~ (SiJiOji ■/hSj2,k) u Z W '* Xk X% (2) 

i,h,h,k 

are space-time points on the lattice.) 
The propagators Sj 2tP = Sj 2y ke WXk can be 
computed from a lattice Dirac equation with a 
momentum source, M being the fermion matrix: 

M it ).Sk,p = e lpXi (momentum source) 

k 

So here the number of matrix inversions is pro- 
portional to the number of momenta. But ev- 
erything else is then there: S p .j 1 — p 7s 
and the quark propagator S(p) — (Sp^u with 
Sp,p = Sfe Sp,kZ lpXk for an amputation and Z q . 

The momentum source method automatically 
performs all the site sums in eq. (2). Another 
possibility, instead of summing over ji (and ji), 
is to choose a particular location for the opera- 
tor, for example setting j2 = 0. Translational 
invariance tells us that this will give the same 
expectation value after averaging over all config- 
urations. For this method we need to solve the 



Dirac equation with a point source at j\ — ji and 
(for extended operators where j± ^ j%) for a small 
number of sources in the region around j% : 

/J Afj.fcfffcjt = h,h (point source) 

k 

(Skj 1 has simply to be fourier transformed to 
S p j 1 then and we have S(p) ~ (Sp.j^jj.) 

For local operators and operators with a small 
number of derivatives the point source method 
would need fewer inversions, but we see from 
Fig.l. that relying on translational invariance to 
carry out the j sums leads to much larger error 
bars. 
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Fi gUX6 1. Z ax i a i with a momentum source method (top 
figure; 20 coring.) and a point source method (120 config.) 

4. RESULTS (WILSON FERMIONS) 

Our calculations are performed on 16 3 32 lat- 
tices at /3 = 6.0 for hopping parameters k — 0.15 15, 
0.1530, 0.1550. If the method works for current 
lattices, one expects to find a high enough /i 2 a 2 
region where the Z factors agree with lattice per- 
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Figure 2 . Z of the scalar and pseudo-scalar operator for 
K = 0.1515, 0.1530, 0.1550 with lower points belonging to 
highest k (lightest quark mass) are shown. 

turbation theory, but are still not destroyed by 
discretisation effects. However the available one- 
loop calculations probably need to be improved 
for g 2 = 6//3 = 1. At low fj, 2 a 2 non-perturbative 
effects and higher order perturbation theory may 
have a big influence as can be seen looking at 
the pseudo-scalar Z g5 (Fig. 2): results for 3 quark 
masses indicate a non-perturbative pion pole con- 
tribution in the chiral limit. The Z factors for 
the operators of the moments of the structure 
functions, e.g. Z( x ) (Fig. 3), seem to show also 
big non-perturbative or higher order perturbation 
theory contributions at low fi 2 a 2 . O(a) discretiza- 
tion errors are in general expected at higher /i 2 a 2 . 
However, for the local vector current (Fig. 3) such 
errors seem to be there at all scales: it is not con- 
stant as expected from its continuum behaviour. 
In the O(a) improved theory it is. 

In a comparison with one-loop lattice pertur- 
bation theory (dashed lines in figures) ||: 



Zo = 1 - g 2 /{l&TT 2 )C F ( l0 lnM + B a ), 
Cf — 4/3 and jo being the anomalous dimen- 
sion, and with tadpole improved theory (solid 
lines in figures) || we find that the improved the- 
ory fits better to the data. For the axial current 
(Fig.l) and the scalar operator (Fig. 2) there is a 
good agreement. We extract Z ax i a i — 0.78. 




Figure 3. Z of the local vector current LVC (re = 0.1530) 
and Z^ (3 re, lower points belong to lowest re) are shown. 
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